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Abstract 

The geometry of a light wavefront, evolving from a initial flat wavefront in the 3-space associated 
with a post-Newtonian relativistic spacetime, is studied numerically by means of the ray tracing 
method. For a discretization of the bidimensional light wavefront, a surface fitting technique is 
used to determine the curvature of this surface. The relationship between the intrinsic curvature of 
the wavefront and the change of the arrival time at different points on the Earth is also numerically 
discussed. 



I. INTRODUCTION 



According to Einstein's theory of general relativity, a generic gravitational field has, at 
a time, three main effects on light: the Shapiro time delay, the bending of rays and the 
curvature of wavefronts. The gravitational field acts respectively as a retarder delaying a 
light wavefront, as a prism tilting a light wavefront, and as a lens curving a light wavefront. 
In special relativity, time delays can be also induced by inertial motion of the observer or the 
source as in the Doppler effect and uniform motion of a observer can result in aberration, 
i.e., an apparent change in the position of the source. 

Hence, the curvature of a initially plane light wavefront by a gravity field is a purely 
general relativistic effect that has no special relativistic analogue. In order to obtain an 
experimental measurement of the curvature of a wave front, Samuel [I] proposed a method 
based on the relation between the differences of arrival time recorded at four points on the 
Earth and the volume of a parallelepiped determined by four points in the curved wavefront 
surface. We will establish a discretized model of the wavefront surface by means of a regular 
triangulation for the study of the curvature(s) of this surface. The main methods and results 
have been recently published in in this work we will only expose a summary of them. 

II. LIGHT PROPAGATION IN A GRAVITATIONAL FIELD 

Let us consider a spacetime {Ai,g) corresponding to a weak gravitational field determined 
by a metric tensor given in a global coordinate system {(z, ct)} by 

gal3 = VcxP + hap, with 7]al3 = diag (1, 1, 1,-1). (l) 

where the coordinate components of the metric deviation haj3 are: 

hab = 2c''^n\\z\\''^6ab, ha4 = -4:C''^K\\z\\'^Za, k^A = ^C'^ K,\\z\\'^ , (2) 

here k := GM represents the gravitational constant of the Sun, located at Z°'{t), and c 
represents the vacuum light speed. 
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The null geodesies, z{t) = (^z{t),t) satisfy the following equations (see [3]): 

Z"" = ^C^h^^^a — [|/i44,t'^fc + hak,t + c{h4a,k — h4k,a)]z'' 
— {h44^l:5f + hak,l — \hkl,a)z^Z^ 

= ga^pz'^z^. 

where, in the first equation, the first and third terms in the right are of order 0(1), while 
the remaining terms are 0(c~^). The second equation is the isotropy constraint satisfied by 
the null geodesies. 

III. NUMERICAL DESCRIPTION OF A SPACELIKE BIDIMENSIONAL LIGHT 
WAVEFRONT 

A. Discretization of the initial wavefront 

We consider a fiat initial surface Sq far from the Sun formed by points (zi, Z2, — C) (with 
C > 0) in an asymptotically Cartesian coordinate system {z}. For the discretization of Sq 
a triangulation is constructed in such a form that each vertex is represented by a complex 
number of the set: 

V := {z = ai + a2'jJ + ascu^ | Oi, 03 G ^, c<j := exp(27ri/3)}, (3) 

The complex plane and the plane iSq may be identified by means of the mapping 
I : z ^ (3fJ(z), Q'(z), — C). Thus a regular triangulation V of the initial wavefront iSq is 
determined. At each vertex in V a photon zq := 2(0) G V with velocity Zq := (0, 0, c) is 
located. The null geodesies equation may be written as a first order differential system 
u = F{u,t), in phase space u = {z,z), which determines in S (the quotient space of A4 
by the global timelike vector field dt associated to the global coordinate system used in the 
post-Newtonian formalism) a flow: z(t) = <ft{zQ, Zq). For each time t the image of Sq under 
the flow $t determines a curved wavefront St- The initial triangulation by V induces a 
triangulation on the final wavefront St, whose vertices we enumerate using the same labels 
used for the corresponding vertices in Sq- 
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Let {y^}^^i be a normal coordinate system with pole at the point P ^ St and associated 
normal reference frame {ejjf^j^. Under the coordinate transformation z'^ y^, from post- 
Newtonian to normal coordinates, the metric tensor 7 on the space S is determined (up to 
terms of first order in F) from 'jab '■= Qab — ^^^f^ by 

^-■=9^5^^"' """""" -'^ = < + A-(i/^-i(ryy/). (4) 

B. Local approximation of the wavefront 

To compute differential magnitudes of the wavefront surface corresponding to the mesh V 
at each inner vertex 2;° G V, j = 1, . . . , J we consider the 1-ring [2;°; 2;jj,]j=i,...,j,fc=o,...,5 formed 
by the six vertices closest to z^: For each 1-ring \z*; z* ] one obtains on the mesh V the 
image 1-ring [zj; Zj^] under the flow ipt- 

In a neighbourhood of the image point zj the wavefront can be approximated by a least- 
squares fitting of the data [zj; zj^,] as the quadric Sf. 

y' = f{y\y') ■■= h^iy')' + (^2yV + hsiy')'. (5) 

using adapted normal coordinates {y^}. 



C. Curvature of the wavefront 

Using coordinates adapted to the quadric St: the metric 7 on £^ induces a metric 

on Sj of the form 

dy^ dy^ , , 

By means of a generalized Gauss formula, the difference of sectional curvatures K and K 
associated with the plane generated by the tangent vectors {vi, V2}, in S and S respectively, 
is the relative sectional curvature i^rei (see [1]): 

i^rel := K{vi, V2) - K{Vi, V2) = A1A2, (7) 

whereas the mean curvature H is determined by half the trace of the second fundamental 
form ir. 

H=l{\, + \2) (8) 
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IV. APPLICATION OF A NUMERICAL INTEGRATOR TO THE CASE OF A 
STATIC SUN 



The ray tracing has been carried out using an integrator based on the classic Taylor 
series method for ordinary differential equations (see [3]). This integrator presents the 
following advantages: allows the control of both the order and the step size employed in 
the method. In this integrator one may use extended arithmetic precision for the highly 
accurate computation required in this problem. (Precision of 120 binary digits and a 
tolerance Tol = l.E — 20 are used to solve this problem.) At each step the null constraint 
equation is preserved. 

Hereafter, we consider the simplest gravitational model generated by a static Sun, con- 
sidered as a point. However, it can be applied to more complex and realistic gravitational 
fields of the solar system, as considered in [6]. We now apply the ray tracing method to a 
light wavefront region propagating along a tubular neighborhood around the 02;^-axis and 
where the inner and outer hexagons have radii of lengths Rq/25 and 2Rq, respectively. 
The initial flat wave surface Sq, perpendicular to the Oz^-axis, is positioned at -100 A.U.. 
The wavefront St is then determined after a trip of 101 A.U., i.e., at the position of the 
Earth. Remember that the focal length for a light ray grazing an opaque Sun is located at 
548.30 A.U. from the Sun. 

In Figure [T| the surface St at the time when the wavefront arrives at the Earth is shown 
using a gray-scale to represent the relative sectional curvature (note we have used a different 
scale on the 02;^-axis). One sees in this figure that the absolute value of the relative 
sectional curvature function defined on St increases as the distance between the photon and 
the 02;^-axis decreases. 

V. VARIATION OF THE TIME OF ARRIVAL AND CURVATURE OF THE 
LIGHT WAVEFRONT 

Suppose four receiving stations located at four points on an Earth hemisphere. The 
arrival time differences between these stations depend on the curvature of the wavefront. 
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FIG. 1: Wavefront surface and relative sectional curvature (gray scale) deformed by a spherical 
gravitational field (a different scale is used for the vertical axis). 



Assumed known the measurements of arrival times corresponding to four points {Qa}a=i 
on the wavefront it is possible to determine an approximation of the wavefront curvature 
in a region far enough from the Sun (say the Earth), without resorting to the ray tracing 
method. 

An estimation of the Gaussian curvature of the wavefront surface can be obtained using 
the notion of the Wald curvature of a metric space established in Distance Geometry (see 
[7]) that in the case of 2-dimensional manifolds agrees with the Gaussian curvature. The 
Wald curvature is determined as the limit of the embedding curvatures of metric quadruples 
isometrically embedded in surfaces of constant curvature (the Euclidean plane M^, the 
2-sphere or the hyperbolic space EP^—^). 

In the hyperbolic plane of curvature — represented by the Blumenthal model 
([?]) we consider the metric quadruple associated with the points {Qa} assuming that the 
geometry of the 3-space in the vicinity of the Earth is Euclidean. 

The Wald curvature associated to the chosen quadruple {Qa} prove to be 

K := = -0.16E-10. (9) 

This result gives an approximation of the total curvature of the wavefront surface under the 
assumption that locally this surface may be identified with a hyperbolic plane in which the 

6 



quadruple considered is isometrically embedded. 
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